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We study the elastic scattering of light in pure KDP crystals and 
KDP crystals with incorporated titanium dioxide nanoparticles. It 
is shown that the optical quality of the crystals decreases insuffi- 
ciently for the used concentrations of nanoparticles. A mathemat- 
ical model of the experimental setup for light scattering measure- 
ments in low-dispersion media is developed and discussed. The 
propagation function of the experimental setup is given in analyt- 
ical form. The relevance of the model is verified with the use of 
experimental scattering data. 



1. Introduction 

The paper is devoted to the optical quahty charac- 
terization of the novel nonlinear optical (NLO) mate- 
rial - potassium dihydrogen phosphate single crystal 
(KH2PO4, KDP) matrix with incorporated titanium 
dioxide nanocrystals (anatase modification) [1]. The hy- 
brid material was designed for the phase matching condi- 
tions and the second harmonic generation efhciency con- 
trol by light due to the resonant excitation of the surface 
states of nanoparticles [2]. We applied the cone-shaped 
interference method to study the effect of Ti02 nanopar- 
ticles on the linear optical response of the medium. A 
special attention was paid to investigate the angular dis- 
tribution of scattered light (scattering indicatrix) and 
scattering extinction losses for a highly transparent low- 
dispersion NLO functional material. 

Experimental data on scattering in low-dispersion ma- 
terials are within a high dynamic range (6-7 orders of the 



signal magnitude). To enable the measurements of weak 
signals of the scattered light, we use a collecting lens 
on the registration CCD array with the angular aper- 
ture that essentially exceeds the initial divergence of a 
laser beam. This leads to an enhancement of the signal- 
to-noise ratio at large scattering angles and simultane- 
ously to the angular averaging of registered data. The 
introduction of the experimental setup apparatus func- 
tion gives possibility to derive the precise angular signal 
distribution using a proper mathematical processing of 
experimental data. 

Previously, a mathematical model was developed in 
[3] . It was shown that the registered scattering signal is 
distorted at small scattering angles. That model can be 
applied to high-dispersive media like that described in 
[4], but it is insufhcient for low-dispersion samples due 
to the lacking of precision. For this reason, a more pre- 
cise model of the experimental setup was developed. The 
consistency of the proposed model is verified, by using 
experimental measurements in series of the single crys- 
tal KDP matrix with different concentrations of Ti02 
nanocrystals [1]. 



2. Experimental Setup 

For isotropic media, the scattering is axially symmetric, 
so it is sufficient to measure a cross-section of the light 
power/intensity angular distribution in any scattering 
plane (Fig. 1). 
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Fig. 1. Scheme of the experimental setup for angular dcijciidencc 
measurements of the scattered light. The lens mounted at a CCD 
array controls the solid angle for the scattering radiation collec- 
tion, 6 denotes the angle between the incidence laser beam and 
scattering signal acquisition directions 

The experimental laser setup for the optical diagnos- 
tics of low-dispersion samples is based on a G-5 goniome- 
ter with different laser sources on its fixed arm: a DPSS 
diode laser (wavelength A = 532 nm. radiation power 
P = 50 mW) or a He-Ne laser (A = 633 nm, P = 10 mW). 
The sample under study is placed at a sample stage, and 
the registration unit is mounted on the movable arm of 
a goniometer at the distance L = 200 mm from the sam- 
ple. It consists of a CCD array AMKO LTI MuLTIray 
(1024 pixels with 25 x 200 /im^ pixel size, 12-bit digital 
resolution) with a lens (diameter d = 9.6 mm) attached 
to it. The rest receiving area of the CCD is completely 
blackened, so that only the light passed through the lens 
reaches the sensitive surface of the CCD array. 

The experimental setup is intended for the investiga- 
tion of a polar angular dependence of the scattered light 
power against to the incidence laser beam propagation 
direction. The magnitudes of radiation power along the 
registration direction 6 in the given solid angle were ob- 
tained via the summation of the all pixels data of the 
CCD array. Due to the focusing lens, the signal on CCD 
is registered on a small number of pixels (^ 30) in each 
measurement, which extremely increases the signal-to- 
noise ratio. 

The indicatrix registered with the use of this experi- 
mental setup is different from the real one due to the spa- 
tial averaging with the lens aperture. So it is reasonable 
to develop a mathematical route to restore the original 
angular distribution from the registered one. The model 
presented here is verified in a direct way, and it can be 
a base for solving the inverse problem in future. 

Scattering data with high angular resolution are re- 
quired for the verification of the model. This was 
achieved due to a reduction of the solid angle for the ra- 
diation collection. The CCD array with a wide aperture 
lens was replaced with a fiber coupled CCD spectrometer 
with the input fiber aperture (d — 410 /im). Unfortu- 
nately, the angular resolution enhancement leads to a 
drastic reduction of the scattering signal level and the 
angular range of data acquisition {6 ~ ±2°). 



The spectroscopic and scattering measurements were 
done with crystalline plates 10 x 10 x 0.8 mm^ cut per- 
pendicular to the optical axis (Z-cut) of the uniaxial 
KDP matrix. A set for the investigation consists of a 
pure KDP single crystalline plate and KDP crystals with 
incorporated Ti02 nanoparticles at different concentra- 
tions (10-^ lO"'', and lO'^ wt.% Ti02) in a mother 
liquor solution. There were two series KDP(Pr) and 
KDP(P) cut from the prismatic {100} and pyramidal 
{101} growth sectors, respectively. The pure KDP crys- 
tals will be denominated simply as P or Pr, and the rest 
of samples will be further designated according to the 
nomenclature in Table 1, where we also give the thick- 
ness h of each sample. In some cases, the Ti02 nanopar- 
ticles concentration is written explicitly in a short form 
(e.g., lO"'^ wt.%), which is referred to the nanoparticles 
concentration in the mother liquor solution. 

Due to the opposite charge signs of the termination 
interfaces of the different growth sectors, the {101} (P) 
growth sector captures more Ti02 nanocrystals in com- 
parison with the {100} (Pr) sector at the growth stage. 
The preliminary research has shown that the average 
concentration of incorporated nanocrystals into the KDP 
matrix is twice less in comparison with the nanoparticles 
concentration in the mother liquor solution. 

3. Mathematical Model, Experimental Results 
and Discussion 

3.1. Optical quality characterization of crystals 

3.1.1. Spectral transmittances 

Optical transmittance spectra of the samples were mea- 
sured with a Perkin Elmer Lambda 35 spectrometer in 
the range 200-1100 nm (the reference is air). The spec- 

Table 1. Spectral properties of the samples (mea- 
sured with a fiber coupled spectrophotometer). "Pr" de- 
notes the prismatic {100} growth sector, "P" — pyramidal 
{101}. Ti02 concentrations are given for mother liquor 
solutions, h is the crystalline plate thickness 



Ti02 cone. 

wt.'/c 


Sample 
iiolaLioii 


/i. 
iiiiii 


Transmittance. % 


A = -hVI mil 


A = ()•->:? um 





P 


0.81 


91.3 ± 1.0 


91.9 ± 0.9 




Pr 


0.83 


92.4 ±0.6 


93.0 ± 0.7 


10-5 


P-5 


0.75 


89.8 ±0.9 


90.2 ± 1.1 




Pr-5 


0.78 


90.3 ± 0.9 


91.0 ± 1.0 


10-4 


P-4 


0.76 


89.8 ± 1.2 


90.8 ± 1.0 




Pr-4 


0.75 


91.4 ± 1.1 


92.0 ± 1.0 


10-3 


Pr-3 


0.81 


89.1 ±0.5 


89.7 ±0.5 
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tra were studied for the samples with pohshed (001) faces 
~0.8 mm in thickness without anireflective coatings. 

The resuhs of these spectral measurements are pre- 
sented in Fig. 2, and the samples notations are given in 
Table 1. The transmittance magnitudes in the visible 
and near-IR ranges are governed, in principle, by the 
reflections from crystalline-air interfaces. 

It is known that that the pure KDP crystals have dif- 
ferent transmittances for prismatic (Pr) and pyramidal 
(P) growth sectors in the UV spectral range [1]. The P 
samples are more transparent due to a less efhcient ab- 
sorption of impurities from the mother liquor. In Fig. 2, 
we show a typical transmittance reduction ^ 75 % at 
A = 270 nm (curve 1). Pr-5 and Pr-3 (curves 2 and 4) 
have transmittances similar to that of the pure one in the 
range 260-350 nm, while Pr-4 (curve 3) demonstrates a 
much higher transmittance in comparison with those of 
the rest ones. It gains ^ 88% at 300 nm. The magni- 
tude is comparable with the transmittance of the pure 
P crystal, and it exceeds the transmission of the pure Pr 
one by about 8%. 

We suggest that the mentioned peculiarity can be ex- 
plained due to the compensation effect of incorporated 
nanoparticles. For the prismatic growth sector, a con- 
centration of 10"'' wt.% Ti02 is not so high for the ef- 
ficient absorption of UV light by anatase nanoparticles 
due to the direct transitions from the valence band to the 
conduction one. However, the concentration of incorpo- 
rated nanoparticles is enough to absorb impurity atoms, 
especially the atoms of polyvalent metals on their devel- 
oped surface during the growth stage. 

We have also tested a homogeneity of the optical trans- 
mittances of crystals over the transverse plane with an 
Ocean Optics type fiber coupled spectrophotometer in 
the range 450-900 nm. The samples were illuminated 
with a broadband coUimated light beam, and the read- 
out of transmittances was performed with a microobjec- 
tive attached to the fiber. It was placed at a distance of 
2 mm from the crystal output interface, so it was possi- 
ble to collect almost all the transmitted photons (both 
ballistic and diffusive ones). We performed 6-10 mea- 
surements for each sample. Their averaged transmission 
coefficients at 532 and 633 nm are presented in Table 1. 

We selected the wavelengths for the characteriza- 
tion of optical homogeneity and scattering proper- 
ties by the following reasons. The available laser 
sources at the mentioned wavelengths correspond to 
different regimes of the excitation of pure and hy- 
brid materials. The wavelength A = 633 nm corre- 
sponds to the nonresonant excitation of the KDP ma- 
trix and Ti02 nanoparticles. The wavelength A = 




200 400 600 800 1000 1200 



Wavelength, nm 

Fig. 2. Spectral transmittances of KDP single crystals with incor- 
porated Ti02 nanoparticles cut from the prismatic growth sector 
(Pr): 1 - pure KDP (Pr), 2- KDP:Ti02 with IQ-^ wt.% (Pr-5), 
3- 10-^ wt.% (Pr-4), 4 - 10-3 wt.% (Pr-3) samples 

532 nm leads to the resonance-type excitation for in- 
trinsic matrix defects [5] and to the electron excita- 
tion from the valence band into the deep level of a de- 
fect in the gap (oxygen vacancies) of anatase nanoparti- 
cles [2]. 

The analysis of the data presented in Table 1 has 
shown that the nonresonant transmittance is a little bit 
higher than that in the resonance excitation case. A 
dispersion of the optical inhomogeneity is the most ev- 
ident for the resonant excitation regime (A = 532 nm) 
and for the lO"'* wt.% Ti02 concentration. This can be 
explained by a high sensitivity of the compensation ef- 
fect to the local concentration of anatase nanoparticles 
in the matrix. 

The comparison of the data in Table 1 with the trans- 
mittance curves in Fig. 1 obtained with two different 
spectroscopic setups shows the high optical transmit- 
tance and the homogeneity of pure KDP and novel hy- 
brid KDP:Ti02 single crystals in the visible and near-IR 
ranges. We observed the 1% difference of the trans- 
mission coefficient for the investigated samples in the vis- 
ible range. In the conventional mode of a Perkin Elmer 
spectrometer (without integrated sphere detection), the 
samples were placed at a distance of 200 mm from the 
1 X 10 mm^ slit of a detector unit. We explain the ob- 
tained minor magnitudes of the Perkin Elmer spectro- 
scopic data by the presence of scattering losses, which 
was proved in the experiment dealing with scattering 
properties of the novel hybrid material. 
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Fig. 3. Scattering indicatrices on the semilogarithmic scale for KDP(P) crystals series at (a) 532-nm and (b) 633-nm wavelengths. 
Squares - freely propagating laser beam; o - pure P crystal; o - KDP:Ti02 P-5 sample; * - KDP:Ti02 P-4 sample 
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Fig. 4. Scattering indicatrices for the series of KDP(P) single 
crystals on the double logarithmic scale at 633 nm. o - pure P 
crystal, o - KDP:Ti02 P-5 sample, * - P-4 sample 

3.1.2. Optical scattering properties 

We have studied the angular distribution of scattered 
hght for both series of KDP(Pr) and KDP(P) crystals in 
two regimes of excitation: (a) resonant regime at 532 nm 
and (b) nonresonant one at 633 nm. The measurements 
were carried out only in the forward hemisphere, i.e. in 
the angle range 0^6*^ 90°. The laser beams propa- 
gated along the optical axis {9 = 0) for Z-axi single crys- 
tals. We performed the axial rotation of the crystalline 
plates and observed the same angular distribution of the 
scattered light (isotropic character of the effect). 

The measured scattering indicatrices for the series of 
KDP(P) crystals are presented at Fig. 3. Squares corre- 
sponds to the angular distribution of the scattered radi- 



ation of the laser beams propagating in air at 532 nm (a) 
and 633 nm (6). The data are presented on the semilog- 
arithmic scale in order to cover a huge dynamic range 
of the described process - about 6-7 orders of the sig- 
nal magnitude reduction at large scattering angles. All 
the data are normalized on the total transmitted power 
Po in the forward hemisphere. It was estimated by the 
incident power of the laser sources and the total trans- 
mittance data presented at Table 1. 

Due to the spatial averaging effect of the lens attached 
to the CCD array, all the curves look like very similar: 
a sharp bell-shaped part and wings. In order to under- 
stand the physical origin of each part, we present the 
data for the KDP(P) series on the double logarithmic 
scale (see Fig. 4). One can clearly see how the presence 
of the lens distorts the shape of the registered indica- 
trix. One can split each curve into three typical ranges, 
where the investigated dependences are distorted at dif- 
ferent extents. 

The first one (0 ^ 1°) is a plateau. In this "blurred" 
interval, the lens collected almost the whole laser beam 
transmitted through the sample - the so-called ballistic 
photon registration geometry. Thus, the acquired data 
do not depend on the angular position of the lens. It is 
the range of the biggest distortions, because, according 
to the spatial averaging, no peculiarities in the indicatrix 
can be resolved at such small angles. 

In the second interval 1° ^ 9 ^ 2°, an abrupt decrease 
of the registered signal is observed, when the lens departs 
from the transmitted laser spot area. This corresponds 
to the direct light to shadow angular range scan with 
large experimental errors (up to 20%) due to the angular 
position uncertainty of the measuring unit. 
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In the third interval ^ 2°, the lens collects only 
the light scattered by a sample. For a medium with- 
out spatially periodic structures, the angular distribu- 
tion of the scattered light is smooth on the logarithmic 
scale. Thus, the spatial averaging of the collecting lens 
has no essential impact on the angular resolution of a 
registered signal. A slight reduction of the angular res- 
olution is compensated by the essential enhancement of 
a registered signal level. 

In the third interval (Fig. 4), we observed an almost 
linear decrease of the registered signal on the double log- 
arithmic scale. This corresponds to the wing parts of the 
same dependences plotted on the semilogarithmic scale 
in Fig. 3. The slope -2.0 ± 0.3 (2° ^ 6* < 45°) is the 
same for the different series of samples and for the differ- 
ent readout wavelengths. This means that the scattered 
radiation power along the 9 direction obeys the inverse 
quadric law AP{9) ~ 0-"^. 

Due to negligible distortions of the angular distribu- 
tion of scattering data in the mentioned intervale, the 
registered indicatrix is very similar to the real one. This 
enables us to estimate the scattering losses in the for- 
ward hemisphere Pscat for all samples under study. The 
method of calculation is based on the integration of the 
scattered light power AP{9) in the solid angle AO over 
the part of the forward hemisphere that corresponds to 
the 3'^'' scattering area in Fig. 4, 

where the lens incoming solid angle Af2 = 
47rsin^ ^lens/2 = 1.8 x 10"^ sr, 6'i„in ~ 1-8° exceeds the 
lens angular aperture ^lens = arctg(d/2L) ~ 1.4° due to 
a finite size of the transmitted beam at the registration 
plane. The data were normalized on the total power Pq 
transmitted into the forward hemisphere. 

The results of estimations of the scattering losses Pscat 
via Eq. (1) are presented in Table 2. The pure P and Pr 
single crystals scatter about 1% of the light radiation for 
both laser wavelengths. The ten-percent increase of the 



Table 2. Scattering losses into the forward hemisphere 
for KDP:Ti02 crystals series normalized on total power 
Pq transmitted into the forward hemisphere, % 





Pure 
KDP 


Ti02 


cone, wt.% 


10-5 


10-* 


10-3 


P (633 nm) 


1.0 


1.7 




1.8 




Pr (633 nm) 


1.1 


2.6 




2.6 


2.6 


P (532 nm) 


1.0 


1.7 




2.5 




Pr (532 nm) 


1.1 


1.7 




1.3 


2.8 



losses in the Pr crystal in comparison with the P one can 
be attributed to a higher concentration of impurities in 
the former [1]. The incorporation of Ti02 nanoparticles 
insignificantly reduces the optical quality of the crystals. 
Even for 10-^ Ti02 wt.% concentration, the scattering 
losses do not exceed 3% that makes it possible to utilize 
the hybrid crystals for the different optical applications. 

In order to explain the obtained data on scattering 
losses, we should refer to the structural characteriza- 
tion of the hybrid media. The high-resolution X-ray 
diffraction analysis has shown that the incorporation of 
nanoparticles has no essential effect on the structural 
perfection of composite materials [1]. In the growth pro- 
cess, the KDP matrix can capture nanoparticles between 
the growth layer stacks, by effectively forming a ID lay- 
ered macroscopic structure . . . KDP:Ti02:KDP. . . with 
a spatial period of 20-30 /xm. This gives possibility to 
incorporate Ti02 nanocrystals, whose nonlinear optical 
response can be controlled by the resonant excitation of 
surface states [2]. 

Thus, the anatase nanocrystals are distributed on the 
planes of the boundaries of KDP growth layer stacks 
with a definite spatial period. The orientations of these 
planes relative to the optical axis are different for the 
P and Pr series. In the Z-cut plates, the optical axis is 
parallel to the planes of boundaries with Ti02 nanopar- 
ticles for the Pr crystals, and it makes an angle of about 
44° for the P one. 

For the nonresonant readout case, the scattering losses 
at 633 nm in both series of KDP(P) and KDP(Pr) crys- 
tals are almost independent of the Ti02 concentration. 
A reduction of the scattering losses by ~ 30% in the 
P series can be explained by the tilt of the planes with 
Ti02 nanoparticles toward the laser beam propagation 
(optical axis) direction, which leads to a decrease of the 
interaction length of light with nanoparticles. 

In the case of the resonant excitation at 532 nm, the 
effect from Ti02 nanoparticles is evident. It is clearly 
seen in Fig. 3 for the P series at the part with indicatrix 
wings. The scattering losses at 532 nm (a) are higher 
for the crystals with incorporated nanoparticles as com- 
pared with those for the pure one. But, in the case of the 
633-nm irradiation (6), they are much alike. The noise 
level reference of the acquisition unit corresponds to the 
wings of the laser beam scattering indicatrix. 

The monotone growth of the scattering losses with the 
concentration of nanoparticles for the P series is deter- 
mined by the following. The preliminary research has 
shown that the KDP:Ti02 single crystals manifest a self- 
focusing effect with the CW laser excitation at 532 nm 
due to the resonant excitation of anatase nanoparticles. 



ISSN 2071-0186. Ukr. J. Phys. 2010. Vol. 55, No. 8 



879 



V.YA. GAYVORONSKY, V.N. STARKOV, M.A. KOPYLOVSKY et al. 




a) h) 

Fig. 5. Cone-shaped interference images for Z-cvA plates of (a) 
pure KDP{P), (6) KDP:Ti02 {P-4) single crystal plates 

The effect leads to the optical contrast enhancement in 
a vicinity of the nanocrystal and to an increase of its 
scattering cross-section. 

The same effect in conjunction with a specific align- 
ment of nanoparticles toward the beam propagation di- 
rection produces a more complicate response to a varia- 
tion of the concentration of nanoparticles in the Pr series. 
An increment of the refractive index at the boundaries of 
KDP growth layer stacks with Ti02 nanoparticles leads 
to the waveguiding effect along the beam propagation 
direction inside KDP layers. The photoinduced waveg- 
uiding effect reduces the scattering of light: at least by 
30% for Pr-5 and by a factor of two for Pr-4 (see Table 2) 
as compared with the nonresonant excitation case. 

The Pr-4 sample possesses unique useful optical prop- 
erties among the rest investigated KDP crystals with in- 
corporated Ti02 nanoparticles: enhanced transmittance 
in the UV-visiblc range and low scattering losses at the 
532-nm resonant laser excitation. 

3.1.3. Cone-shaped interference analysis 

Pure KDP belongs to uniaxial negative crystals. The 
incorporation of Ti02 nanoparticles produces additional 
stresses/strains in crystals. Due to this fact, the crys- 
tals acquire the anomalous biaxiality (AB) which can 
be measured for Z-cut plates within the cone-shaped in- 
terference analysis. The technique is the following: a 
sample is mounted between crossed polarizers, and it is 
illuminated with a cone-shaped beam. 

The cone-shaped interference images for Z-cut plates 
of pure KDP(P) (a) and P-4 (6) samples are presented 
in Fig. 5. For isotropic media like pure KDP(P), two 
isogyres overlap along the optical axis and form the dark 
cross (a). The incorporation of nanoparticles induces 
the AB effect and a lateral optical anisotropy. The last 
one exhibits itself as the isogyre divergence (6). The 



extent of the AB effect is determined from the magnitude 
of isogyre divergence (the angle 2V) in a cone-shaped 
interference pattern. 

The distance between the isogyres is independent of 
the rotational orientation of a sample and depends on 
the 2V value of a material. When this 2V value is large 
(about 40 50°), two isogyres are rarely seen at a single 
image. However, for materials with a weak AB effect, 
the 2V value is small and can be measured directly from 
the image. 

Values of 2V were measured for KDP(P) crystals over 
different lateral points of samples. The results are sum- 
marized in Table 3. They show that the strains in pure 
KDP single crystals are very small. The manifestation 
of the AB effect enhances with increase in the concentra- 
tion of incorporated Ti02 nanoparticles in the matrix. 
However, even for 10"^ Ti02 wt.%, the value of 2V is 
less than 20', which indicates a good optical quality of 
the samples under study. 

3.2. The mathematical model 

As was shown above, the registered scattering indicatri- 
ces are distorted by the presence of the lens, especially 
in areas 1 and 2 in Fig. 4. Let 6*0 denote the current 
position of the lens center. Then the registered depen- 
dences P{0o) are quite precise at > 2°. But, at small 
angles, they need an additional processing. When the 
"wings" of indicatrices are compared on the logarithmic 
scale, the central part |^o| ^ 2° should be given on a 
linear scale for convenience. 

There are two possible ways to obtain the real view of 
the indicatrix at ^ 2°. The first way is experimental 
(using an optical fiber spectrophotometer instead of a 
CCD with a lens), and the second one is a way of math- 
ematical processing. The first way is rather complicated 
and requires a lot of time to execute experiments. So it 
is reasonable to develop a simulation procedure to solve 
this problem, and the current mathematical model is the 
first step on this way. The verification of the model was 
performed, by using a 532-nm continuous laser and the 
series of KDP(P) crystals. 



Table 3. Results of the cone-shaped interference 
analysis for KDP(P) series 



Ti02 cone, 
wt.% 


Intrinsic strain anomalies 

(number of points) 


2yav. 


0-56' 


5-10' 


10-15' 


15-20' 





60 


40 






6 


10-5 






80 


20 


12 


10-4 






40 


60 


17 
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Fig. 6. Lines of equal radiation intensities on the lens: when the laser beam axis is in (a) and beyond (6) the lens aperture 



For isotropic samples, the intensity distribution of 
scattered radiation is axially symmetric, so the lines with 
the same intensity are circumferences of radii p = L sin 6 
(Fig. 6), where 9 is the azimuth angle, and L is the dis- 
tance between the sample and the lens (Fig. 1). As the 
registration unit moves along a circle, it is reasonable to 
use a spherical coordinate system. But, to process the 
data registered at the rotation angles |6'o| ^ 2° of a mov- 
able goniometer arm, a cylindric coordinate system can 
be used. The coordinate origin O is the point where the 
laser beam axis hits the center of a lens at 6q = 0. Then, 
at 7^ 0, the lens center coordinate x = Lsin^o ~ LOq 
(Fig. 6), a is the polar angle, and r is the distance from 
an arbitrary point on the lens aperture to the lens cen- 
ter. The laser beam axis in Fig. 6 is perpendicular to 
the plane of the figure. 

For convenience, u{.r) in the mathematical part of the 
work means the registered scattered power AP{9o), and 
i'{p) stands for the real scattered intensity /scat(^)- Then 
the infinitesimal amount of a registered radiation power 
at each point of the lens can be given by 

Au^i^{p)G{x,p,a)AS, (2) 

where G is the propagation function at a specific point on 
the lens, and AS is an area element of the lens aperture 
(Fig. 6). 

The integration of Eq. (2) over the lens aperture gives 
the registered radiation power at a certain position 6q of 
a CCD with a lens. If we assume that the propagation 
function depends only on the polar radius of a point on 
the lens G = G{r), we can write 

P2{x) a2{x,p) 

u{x) = J pi'{p) j G{x, p,a)dpda, x£[—a,a]. (3) 

Here, G{x, p, a) = G{r) according to 

= p^ + x^ — 2px cos a. (4) 



According to the geometry of the experiment (L = 
200 mm, a ^ 7 mm), Eqs. (3)-(4) are written with an 
error less than 0.08 %. 

The dependence G{r) can be directly measured if an 
ideally thin laser beam scans the lens along its diame- 
ter. In this case, at a point r = Rq {Rq = d/2 is the 
lens radius), the abrupt drop in a registered signal from 
a finite value to zero would be observed. But, by math- 
ematical reasons, G{r) should be continuous at every 
point of r G [0,-|-oo), and this abrupt drop should be 
eliminated. 

As the laser beam diameter 5 « 0.9 mm is sufficiently 
less than the lens size, the experimental data for a freely 
propagating laser beam were taken as a propagation 
function. It can be approximated with a continuous 
function, and its form is similar to the real propagation 
function because 5 <^d. The results of the experimental 
measurements are shown in Fig. 7. As a mathematical 
approximation, we use the function 

G{r)=x{Rl-r^){AeM{r/rAf)- 



-{rlTBf - Cexp(-Y3^^^)x(ri - r^)), (5) 
where i?i = 5.509 mm, 

A = 2.223, TA = 5.808 mm, G = 1.277, 

rs = 2.357 mm, rc = 3.574 mm, tr = 4.711 mm, 

and x(s) is the Heaviside function equal to unity for 

s ^ and equal to zero for ,s < 0. 

In Eq. (5), we have G(r) = for r > Ri. Here, 
Ri > Ro, because the beam radius is non-zero. The 
relative approximation error is 3.6%, while the relative 
experimental error is « 5%. 
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Fig. 7. Propagation function of the experimental setup (solid line) 
and experimental data for a freely propagated beam (dots) 

From the symmetry in Fig. 6, it is clear that, in 
Eq. (3), ai{x,p) = —a2{x,p). Using Eq. (4), we can 
write 



arccos 



V 2pa; ) 



p + X > Ri, 
p + X ^ Rl. 



(6) 



Using Eq. (6) and the function G(x, p, a), one can found 
the inner integral in Eq. (3). As G{x, p, a) is a symmetric 
function, it can be written as 



"2(2;, p) 



«2(x,p) 



K{x,p) = ^ J G{x,p,a)da= J G{x,p,a)da. (7) 

ai{x,p) 

Equation (7) enables us to write the integral relationship 
Eq. (3) in the form 



i{x) = 2 J pu{p)K{x,p)dp, xG[—a,a]. 



(8) 



pi (a:) 



As for determining the integration limits pi{x) and 
P2{x), two positions of the lens with respect to the cen- 
ter of a laser beam are possible. In the first case, the 
center of the beam lies inside the circle of the radius Ri 
(Fig. 6, a, X ^ Rl). In the second case, the center of the 
beam is outside the circle (Fig. 6,&, a; > At first, 
we discuss both of the cases apart. 

1) a; ^ Rl. 

According to Eq. (6), in this case, the lines of equal 
radiation intensities inside the circle with Ri may be of 
two types: 



a) Circumferences of radii ^ p ^ (i?i — x). Then 
a2{x, p) = TT and 



7r 

Ko{x,p) = j G{x,p,a)da, 


Ri—x 

u{x) = 2 j pL'{p)Ko{x,p)dp. 



(9) 



b) Arcs of the circumferences {Ri — x) < p < {Ri + x). 
Then 0:2(0;,^) < tt according to Eq. (6) and 

0!2(a;,p) 



Ki{x,p) = j G{x, p,a)da, 



Ri+x 

u{x) = 2 j pv{p)Ki{x, p)dp. 

Ri—x 



(10) 



The functions Kq{x,p) and Ki{x,p) arc presented in 
Fig. 8. They can be united into a 2D continuous function 
Kix,p). 

2) x> Rl. 

In this case, the lines of equal radiation intensities are 
also arcs of the circumferences with radii {x — Ri) < p < 
(x + Ri). Here, u{x) can be given identically to Eq. (10), 
but with different integration limits. The total integra- 
tion relation Eq. (8) can be written using Eqs. (9)-(10) 
and Eq. (5) as follows: 

2Ri 

u{x) = 2 J pi'{p)K{x, p)dp, xG[—a,a\, (11) 



where K{x, p) is the union of the functions Kq{x, p) and 
Ki{x,p) (Fig. 8). 

There are two ways to interpret Eq. (11). When the 
real intensity distribution of the scattered light i'{p) is 
known, Eq. (11) becomes the integration forrrmla one can 
use to mathematically obtain the resulting data of the 
experiment with certain parameters of the laser setup. 
In this case, the direct problem is solved, and no com- 
putational difficulties arise. 

Another situation appears at solving the inverse prob- 
lem. Here, u{x) is known, and the task is to recon- 
struct the real intensity distribution v{p) of the radi- 
ation scattered by the sample. Then Eq. (11) serves 
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K[x,p] 




Fig. 8. Form of the integrand function K{x,p) 

as the Fredholm first-kind integral equation, and its 
solution is a solution of an incorrect problem. In 
this case, special regularization methods are required 
[6, 7]. 

At the stage of verification of the consistency of the 
mathematical model of the experimental laser setup for 
the optical diagnostics of low-dispersion samples, it is 
suitable to have the experimental data for ^{p) and u{x). 
For this purpose, we studied 

a) Freely propagated 532-nm laser beam; 

b) Laser beam passes through a pure KDP(P) crystal; 

c) Laser beam passes through a KDP crystal doped 
with Ti02 nanoparticles (sample P-4). 

To experinicritally obtain the data on the real inten- 
sity distribution we used the experimental setup 
with an optical-fiber (diameter d = 410 fim) spectropho- 
tometer instead of a CCD with a lens. Results of the 
experiment reveal that, in all three investigated cases, 
the propagating laser beam has the Gaussian intensity 
distribution 

u{p)=eM-ip/pof)- (12) 

A beam size reduction in the registration plane after 
the transmission of samples at 532 nm was shown (po = 
0.447 mm in the case of free propagation of the beam) . In 
case of KDP with Ti02 (0.421 mm), it is less than that 
for pure KDP (0.432 mm). This can be explained by 
the photoinduced self-focusing effect in these samples at 
the 532-nm CW laser excitation, which will be published 
elsewhere in detail. This result proves that knowing the 
real scattering indicatrix can be essential for the optical 
characterization of a sample. 



u(x)-10^ 




X, mm 



Fig. 9. Comparison of the experimental data at 532 nm and the 
simulation results, (a) Dots — experimental data for P-4 sam- 
ple, solid line - simulation, dashed line — precise measurements 
(normalized for convenience). (6) Simulation results: 1 - freely 
propagated laser beam, 2 - KDP{P) pure, 3 - P-4 sample 

In Fig. 9, a, the experimental data (dots) and the cal- 
culated u{x) dependence (solid line) for P-4 sample are 
presented. We also show the precise scattering data for 
P-4 sample (dashed line, normalized to ~ 5 x 10~^ for 
convenience). Figure 9,6 shows the difference between 
the calculated indicatrices for a freely propagated laser 
beam, pure KDP(P), and P-4 sample. The result pre- 
sented is in good agreement with spectral data (Table 1). 
The dependences u{x) for all the samples are obtained, 
by using Eq. (11) and functions (12), the approximation 
error being less than 5%. 

Based on the results of the work, it can be stated that 
the mathematical model of the experimental laser setup 
for the optical diagnostics of low-dispersion samples is 
consistent, and it can serve as a base for developing a 
way to solve the inverse problem, i.e. to get the real 
scattering data /scat (6*) out of the registered AP(^o) us- 
ing mathematical calculations. 

4. Summary 

In this work, the optical characterization of series of 
KDP crystals doped with Ti02 nanoparticles (anatase 
modification) in various concentrations was performed. 
The characterization includes measurements of scatter- 
ing indicatrices and optical transmission spectra and the 
cone-shaped interference method. 

The investigations show that the incorporation of 
Ti02 nanoparticles slightly changes the optical quality 
of KDP crystals. All of the characterization results con- 
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tribute to a good optical quality of the investigated sam- 
ples: transmittances in the visible range ~ 90%, scat- 
tering losses < 3%, and the anomalous biaxiality values 
2V < 20'. 

A mathematical model of the experimental setup was 
developed in order to enhance the angular resolution in 
the scattering indicatrices at angles 0^2°, where the 
registered data arc blurred due to the presence of a lens 
in the measuring circuit. The model was verified in the 
direct way and showed itself to be consistent with an 
error less than 5%. 

The work was supported by STCU-NASU grant 
No. 4956. 
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